Profiling the Interaction between Human Serum Albumin and Clinically Relevant HIV Reverse Transcriptase Inhibitors.

Tenofovir (TFV) is the active form of the prodrugs tenofovir disoproxil fumarate (TDF) and tenofovir alafenamide (TAF), both clinically prescribed as HIV reverse transcriptase inhibitors. The biophysical interactions between these compounds and human serum albumin (HSA), the primary carrier of exogenous compounds in the human bloodstream, have not yet been thoroughly characterized. Thus, the present study reports the interaction profile between HSA and TFV, TDF, and TAF via UV-Vis, steady-state, and time-resolved fluorescence techniques combined with isothermal titration calorimetry (ITC) and in silico calculations. A spontaneous interaction in the ground state, which does not perturb the microenvironment close to the Trp-214 residue, is classified as weak. In the case of HSA/TFV and HSA/TDF, the binding is both enthalpically and entropically driven, while for HSA/TAF, the binding is only entropically dominated. The binding constant (Ka) and thermodynamic parameters obtained via ITC assays agree with those obtained using steady-state fluorescence quenching measurements, reinforcing the reliability of the data. The small internal cavity known as site I is probably the main binding pocket for TFV due to the low steric volume of the drug. In contrast, most external sites (II and III) can better accommodate TAF due to the high steric volume of this prodrug. The cross-docking approach corroborated experimental drug-displacement assays, indicating that the binding affinity of TFV and TAF might be impacted by the presence of different compounds bound to albumin. Overall, the weak binding capacity of albumin to TFV, TDF, and TAF is one of the main factors for the low residence time of these antiretrovirals in the human bloodstream; however, positive cooperativity for TAF and TDF was detected in the presence of some drugs, which might improve their residence time (pharmacokinetic profile).

Synthesis and anti-HIV activities of phorbol derivatives.

In this study, 37 derivatives of phorbol esters were synthesized and their anti-HIV-1 activities evaluated, building upon our previous synthesis of 51 phorbol derivatives. 12-Para-electron-acceptor-trans-cinnamoyl-13-decanoyl phorbol derivatives stood out, demonstrating remarkable anti-HIV-1 activities and inhibitory effects on syncytia formation. These derivatives exhibited a higher safety index compared with the positive control drug. Among them, 12-(trans-4-fluorocinnamoyl)-13-decanoyl phorbol, designated as compound 3c, exhibited the most potent anti-HIV-1 activity (EC50 2.9 nmol·L-1, CC50/EC50 11 117.24) and significantly inhibited the formation of syncytium (EC50 7.0 nmol·L-1, CC50/EC50 4891.43). Moreover, compound 3c is hypothesized to act both as an HIV-1 entry inhibitor and as an HIV-1 reverse transcriptase inhibitor. Isothermal titration calorimetry and molecular docking studies indicated that compound 3c may also function as a natural activator of protein kinase C (PKC). Therefore, compound 3c emerges as a potential candidate for developing new anti-HIV drugs.

Targeting the Structural Maturation Pathway of HIV-1 Reverse Transcriptase.

Formation of active HIV-1 reverse transcriptase (RT) proceeds via a structural maturation process that involves subdomain rearrangements and formation of an asymmetric p66/p66' homodimer. These studies were undertaken to evaluate whether the information about this maturation process can be used to identify small molecule ligands that retard or interfere with the steps involved. We utilized the isolated polymerase domain, p51, rather than p66, since the initial subdomain rearrangements are largely limited to this domain. Target sites at subdomain interfaces were identified and computational analysis used to obtain an initial set of ligands for screening. Chromatographic evaluations of the p51 homodimer/monomer ratio support the feasibility of this approach. Ligands that bind near the interfaces and a ligand that binds directly to a region of the fingers subdomain involved in subunit interface formation were identified, and the interactions were further characterized by NMR spectroscopy and X-ray crystallography. Although these ligands were found to reduce dimer formation, further efforts will be required to obtain ligands with higher binding affinity. In contrast with previous ligand identification studies performed on the RT heterodimer, subunit interface surfaces are solvent-accessible in the p51 and p66 monomers, making these constructs preferable for identification of ligands that directly interfere with dimerization.

Exploring novel HIV-1 reverse transcriptase inhibitors with drug-resistant mutants: A double mutant surprise.

HIV-1 reverse transcriptase (RT) remains a key target for HIV drug development. As successful management of the disease requires lifelong treatment, the emergence of resistance mutations is inevitable, making development of new RT inhibitors, which remain effective against resistant variants crucial. To this end, previous computationally guided drug design efforts have resulted in catechol diether compounds, which inhibit wildtype RT with picomolar affinities and appear to be promising preclinical candidates. To confirm that these compounds remain potent against Y181C, a widespread mutation conferring resistance to first generation inhibitors, they were screened against the HIV-1 N119 clinical isolate, reported as a Y181C single mutant. In comparison to a molecular clone with the same mutation, N119 appears less susceptible to inhibition by our preclinical candidate compounds. A more detailed sequencing effort determined that N119 was misidentified and carries V106A in combination with Y181C. While both indolizine and naphthalene substituted catechol diethers are potent against the classical Y181C single mutant, the addition of V106A confers more resistance against the indolizine derivatives than the naphthalene derivatives. Crystal structures presented in this study highlight key features of the naphthyl group, which allow these compounds to remain potent in the double mutant, including stronger interactions with F227 and less reliance on V106 for stabilization of the ethoxy-uracil ring, which makes critical hydrogen bonds with other residues in the binding pocket.

Strategies in the Design and Development of Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs).

AIDS (acquired immunodeficiency syndrome) is a potentially life-threatening infectious disease caused by human immunodeficiency virus (HIV). To date, thousands of people have lost their lives annually due to HIV infection, and it continues to be a big public health issue globally. Since the discovery of the first drug, Zidovudine (AZT), a nucleoside reverse transcriptase inhibitor (NRTI), to date, 30 drugs have been approved by the FDA, primarily targeting reverse transcriptase, integrase, and/or protease enzymes. The majority of these drugs target the catalytic and allosteric sites of the HIV enzyme reverse transcriptase. Compared to the NRTI family of drugs, the diverse chemical class of non-nucleoside reverse transcriptase inhibitors (NNRTIs) has special anti-HIV activity with high specificity and low toxicity. However, current clinical usage of NRTI and NNRTI drugs has limited therapeutic value due to their adverse drug reactions and the emergence of multidrug-resistant (MDR) strains. To overcome drug resistance and efficacy issues, combination therapy is widely prescribed for HIV patients. Combination antiretroviral therapy (cART) includes more than one antiretroviral agent targeting two or more enzymes in the life cycle of the virus. Medicinal chemistry researchers apply different optimization strategies including structure- and fragment-based drug design, prodrug approach, scaffold hopping, molecular/fragment hybridization, bioisosterism, high-throughput screening, covalent-binding, targeting highly hydrophobic channel, targeting dual site, and multi-target-directed ligand to identify and develop novel NNRTIs with high antiviral activity against wild-type (WT) and mutant strains. The formulation experts design various delivery systems with single or combination therapies and long-acting regimens of NNRTIs to improve pharmacokinetic profiles and provide sustained therapeutic effects.

Human immunodeficiency virus 1 5'-leader mutations in plasma viruses before and after the development of reverse transcriptase inhibitor-resistance mutations.

Human immunodeficiency virus 1 (HIV-1) reverse transcriptase (RT) initiation depends on interaction between viral 5'-leader RNA, RT and host tRNA3Lys. Therefore, we sought to identify co-evolutionary changes between the 5'-leader and RT in viruses developing RT-inhibitor resistance mutations. We sequenced 5'-leader positions 37-356 of paired plasma virus samples from 29 individuals developing the nucleoside RT inhibitor (NRTI)-resistance mutation M184V, 19 developing a non-nucleoside RT inhibitor (NNRTI)-resistance mutation and 32 untreated controls. 5'-Leader variants were defined as positions where ≥20 % of next-generation sequencing (NGS) reads differed from the HXB2 sequence. Emergent mutations were defined as nucleotides undergoing a ≥4-fold change in proportion between baseline and follow-up. Mixtures were defined as positions containing ≥2 nucleotides each present in ≥20 % of NGS reads. Among 80 baseline sequences, 87 positions (27.2 %) contained a variant; 52 contained a mixture. Position 201 was the only position more likely to develop a mutation in the M184V (9/29 vs 0/32; P=0.0006) or NNRTI-resistance (4/19 vs 0/32; P=0.02; Fisher's exact test) groups than the control group. Mixtures at positions 200 and 201 occurred in 45.0 and 28.8 %, respectively, of baseline samples. Because of the high proportion of mixtures at these positions, we analysed 5'-leader mixture frequencies in two additional datasets: five publications reporting 294 dideoxyterminator clonal GenBank sequences from 42 individuals and six National Center for Biotechnology Information (NCBI) BioProjects reporting NGS datasets from 295 individuals. These analyses demonstrated position 200 and 201 mixtures at proportions similar to those in our samples and at frequencies several times higher than at all other 5'-leader positions. Although we did not convincingly document co-evolutionary changes between RT and 5'-leader sequences, we identified a novel phenomenon, wherein positions 200 and 201 immediately downstream of the HIV-1 primer binding site exhibited an extraordinarily high likelihood of containing a nucleotide mixture. Possible explanations for the high mixture rates are that these positions are particularly error-prone or provide a viral fitness advantage.

Charge Engineering of the Nucleic Acid Binding Cleft of a Thermostable HIV-1 Reverse Transcriptase Reveals Key Interactions and a Novel Mechanism of RNase H Inactivation.

Coupled with PCR, reverse transcriptases (RTs) have been widely used for RNA detection and gene expression analysis. Increased thermostability and nucleic acid binding affinity are desirable RT properties to improve yields and sensitivity of these applications. The effects of amino acid substitutions in the RT RNase H domain were tested in an engineered HIV-1 group O RT, containing mutations K358R/A359G/S360A and devoid of RNase H activity due to the presence of E478Q (O3MQ RT). Twenty mutant RTs with Lys or Arg at positions interacting with the template-primer (i.e., at positions 473-477, 499-502 and 505) were obtained and characterized. Most of them produced significant amounts of cDNA at 37, 50 and 65 °C, as determined in RT-PCR reactions. However, a big loss of activity was observed with mutants A477K/R, S499K/R, V502K/R and Y505K/R, particularly at 65 °C. Binding affinity experiments confirmed that residues 477, 502 and 505 were less tolerant to mutations. Amino acid substitutions Q500K and Q500R produced a slight increase of cDNA synthesis efficiency at 50 and 65 °C, without altering the KD for model DNA/DNA and RNA/DNA heteroduplexes. Interestingly, molecular dynamics simulations predicted that those mutations inactivate the RNase H activity by altering the geometry of the catalytic site. Proof of this unexpected effect was obtained after introducing Q500K or Q500R in the wild-type HIV-1BH10 RT and mutant K358R/A359G/S360A RT. Our results reveal a novel mechanism of RNase H inactivation that preserves RT DNA binding and polymerization efficiency without substituting RNase H active site residues.

5-Nitro-3-(2-(4-phenylthiazol-2-yl)hydrazineylidene)indolin-2-one derivatives inhibit HIV-1 replication by a multitarget mechanism of action.

In the effort to identify and develop new HIV-1 inhibitors endowed with innovative mechanisms, we focused our attention on the possibility to target more than one viral encoded enzymatic function with a single molecule. In this respect, we have previously identified by virtual screening a new indolinone-based scaffold for dual allosteric inhibitors targeting both reverse transcriptase-associated functions: polymerase and RNase H. Pursuing with the structural optimization of these dual inhibitors, we synthesized a series of 35 new 3-[2-(4-aryl-1,3-thiazol-2-ylidene)hydrazin-1-ylidene]1-indol-2-one and 3-[3-methyl-4-arylthiazol-2-ylidene)hydrazine-1-ylidene)indolin-2-one derivatives, which maintain their dual inhibitory activity in the low micromolar range. Interestingly, compounds 1a, 3a, 10a, and 9b are able to block HIV-1 replication with EC50 < 20 µM. Mechanism of action studies showed that such compounds could block HIV-1 integrase. In particular, compound 10a is the most promising for further multitarget compound development.

Alkaloids as Potential Anti-HIV Agents.

BACKGROUND: Alkaloids are nitrogen-containing compounds that are naturally occurring and have a variety of biological activities, including antimicrobial properties. In this study, the authors used a molecular docking approach to evaluate the anti-HIV potential of 64 alkaloids. METHODS: The authors used the Molegro Virtual Docker software to dock the alkaloids into the active sites of three HIV enzymes: protease, integrase, and non-nucleoside reverse transcriptase (NNRT). The docking scores were used to assess the potential of the alkaloids to inhibit the enzymes. RESULTS: The results showed the alkaloids to have good potential to inhibit the enzymes. Tubocurarine and reserpine were found to be the most potent alkaloids, with docking scores of -123.776 and - 114.956, respectively. CONCLUSION: The authors concluded that tubocurarine and reserpine could be further promoted as potential lead molecules for the development of new anti-HIV drugs.

Cryo-EM structures of wild-type and E138K/M184I mutant HIV-1 RT/DNA complexed with inhibitors doravirine and rilpivirine.

Structures trapping a variety of functional and conformational states of HIV-1 reverse transcriptase (RT) have been determined by X-ray crystallography. These structures have played important roles in explaining the mechanisms of catalysis, inhibition, and drug resistance and in driving drug design. However, structures of several desired complexes of RT could not be obtained even after many crystallization or crystal soaking experiments. The ternary complexes of doravirine and rilpivirine with RT/DNA are such examples. Structural study of HIV-1 RT by single-particle cryo-electron microscopy (cryo-EM) has been challenging due to the enzyme's relatively smaller size and higher flexibility. We optimized a protocol for rapid structure determination of RT complexes by cryo-EM and determined six structures of wild-type and E138K/M184I mutant RT/DNA in complexes with the nonnucleoside inhibitors rilpivirine, doravirine, and nevirapine. RT/DNA/rilpivirine and RT/DNA/doravirine complexes have structural differences between them and differ from the typical conformation of nonnucleoside RT inhibitor (NNRTI)-bound RT/double-stranded DNA (dsDNA), RT/RNA-DNA, and RT/dsRNA complexes; the primer grip in RT/DNA/doravirine and the YMDD motif in RT/DNA/rilpivirine have large shifts. The DNA primer 3'-end in the doravirine-bound structure is positioned at the active site, but the complex is in a nonproductive state. In the mutant RT/DNA/rilpivirine structure, I184 is stacked with the DNA such that their relative positioning can influence rilpivirine in the pocket. Simultaneously, E138K mutation opens the NNRTI-binding pocket entrance, potentially contributing to a faster rate of rilpivirine dissociation by E138K/M184I mutant RT, as reported by an earlier kinetic study. These structural differences have implications for understanding molecular mechanisms of drug resistance and for drug design.

Human endogenous retrovirus-K (HERV-K) reverse transcriptase (RT) structure and biochemistry reveals remarkable similarities to HIV-1 RT and opportunities for HERV-K-specific inhibition.

Human endogenous retroviruses (HERVs) comprise nearly 8% of the human genome and are derived from ancient integrations of retroviruses into the germline. The biology of HERVs is poorly defined, but there is accumulating evidence supporting pathological roles in diverse diseases, such as cancer, autoimmune, and neurodegenerative diseases. Functional proteins are produced by HERV-encoded genes, including reverse transcriptases (RTs), which could be a contributor to the pathology attributed to aberrant HERV-K expression. To facilitate the discovery and development of HERV-K RT potent and selective inhibitors, we expressed active HERV-K RT and determined the crystal structure of a ternary complex of this enzyme with a double-stranded DNA substrate. We demonstrate a range of RT inhibition with antiretroviral nucleotide analogs, while classic nonnucleoside analogs do not inhibit HERV-K RT. Detailed comparisons of HERV-K RT with other known RTs demonstrate similarities to diverse RT families and a striking similarity to the HIV-1 RT asymmetric heterodimer. Our analysis further reveals opportunities for selective HERV-K RT inhibition.

In silico design, synthesis and anti-HIV activity of quinoline derivatives as non-nucleoside reverse transcriptase inhibitors (NNRTIs).

A series of quinoline derivatives has been designed, synthesized and screened for their anti-HIV properties. The drug-like properties of compounds were evaluated first and then molecular docking using DS v20.1.0.19295 software showed that the compounds behaved as non-nucleoside reverse transcriptase inhibitors (NNRTIs) while interacting at the allosteric site of target HIV-RT protein (PDB:3MEC). The docking results revealed that all compounds formed hydrogen bonds with Lys101, Lys103, Val179, Tyr188, Gln190, Gly190, Pro225, Phe227, and Tyr318, and showed π-interaction with Tyr188 and Tyr318. TOPKAT (Toxicity Prediction by Komputer Assisted Technology) results confirmed that the compounds were found to be less toxic than the reference drugs. Density functional theory (DFT) analysis was performed to assess the binding affinity of all compounds. Further, molecular dynamics (MD) simulations were performed on compound 6 and delavirdine with HIV-RT enzyme. Comprehensive MD analyses showed a similar pattern of conformational stability and flexibility in both the complexes suggesting alike inhibitory action. The hydrogen-bonding interactions and the binding energy of active-site residues for the compound 6 complex revealed strong inhibitory activity than the reference (delavirdine) complex. Thus, the compound 6 might act as a potential inhibitor against HIV-RT. Overall, this study revealed that compound 6 (5-hydroxy-N-(4-methyl-2-oxo-1,2-dihydroquinolin-8-yl) thiophene-2-sulfonamide) has prudent anti-HIV activity against both HIV-1 (SI = 2.65) and HIV-2 (SI = 2.32) that can further be utilised in drug discovery against HIV virus.

Relative domain orientation of the L289K HIV-1 reverse transcriptase monomer.

HIV-1 reverse transcriptase (RT) is a heterodimer comprised p66 and p51 subunits (p66/p51). Several single amino acid substitutions in RT, including L289K, decrease p66/p51 dimer affinity, and reduce enzymatic functioning. Here, small-angle X-ray scattering (SAXS) with proton paramagnetic relaxation enhancement (PRE), 19 F site-specific NMR, and size exclusion chromatography (SEC) were performed for the p66 monomer with the L289K mutation, p66L289K . NMR and SAXS experiments clearly elucidated that the thumb and RNH domains in the monomer do not rigidly interact with each other but are spatially close to the RNH domain. Based on this structural model of the monomer, p66L289K and p51 were predicted to form a heterodimer while p66 and p51L289K not. We tested this hypothesis by SEC analysis of p66 and p51 containing L289K in different combinations and clearly demonstrated that L289K substitution in the p51 subunit, but not in the p66 subunit, reduces p66/p51 formation. Based on the derived monomer model and the importance of the inter-subunit RNH-thumb domain interaction in p66/p51, validated by SEC, the mechanism of p66 homodimer formation was discussed.

HIV-1 drug resistance profiling using amino acid sequence space cartography.

MOTIVATION: Human immunodeficiency virus (HIV) drug resistance is a global healthcare issue. The emergence of drug resistance influenced the efficacy of treatment regimens, thus stressing the importance of treatment adaptation. Computational methods predicting the drug resistance profile from genomic data of HIV isolates are advantageous for monitoring drug resistance in patients. However, existing computational methods for drug resistance prediction are either not suitable for emerging HIV strains with complex mutational patterns or lack interpretability, which is of paramount importance in clinical practice. The approach reported here overcomes these limitations and combines high accuracy of predictions and interpretability of the models. RESULTS: In this work, a new methodology based on generative topographic mapping (GTM) for biological sequence space representation and quantitative genotype-phenotype relationships prediction purposes was introduced. The GTM-based resistance landscapes allowed us to predict the resistance of HIV strains based on sequencing and drug resistance data for three viral proteins [integrase (IN), protease (PR) and reverse transcriptase (RT)] from Stanford HIV drug resistance database. The average balanced accuracy for PR inhibitors was 0.89 ± 0.01, for IN inhibitors 0.85 ± 0.01, for non-nucleoside RT inhibitors 0.73 ± 0.01 and for nucleoside RT inhibitors 0.84 ± 0.01. We have demonstrated in several case studies that GTM-based resistance landscapes are useful for visualization and analysis of sequence space as well as for treatment optimization purposes. Here, GTMs were applied for the in-depth analysis of the relationships between mutation pattern and drug resistance using mutation landscapes. This allowed us to predict retrospectively the importance of the presence of particular mutations (e.g. V32I, L10F and L33F in HIV PR) for the resistance development. This study highlights some perspectives of GTM applications in clinical informatics and particularly in the field of sequence space exploration. AVAILABILITY AND IMPLEMENTATION: https://github.com/karinapikalyova/ISIDASeq. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structural Basis of 2-Phenylamino-4-phenoxyquinoline Derivatives as Potent HIV-1 Non-Nucleoside Reverse Transcriptase Inhibitors.

New target molecules, namely, 2-phenylamino-4-phenoxyquinoline derivatives, were designed using a molecular hybridization approach, which was accomplished by fusing the pharmacophore structures of three currently available drugs: nevirapine, efavirenz, and rilpivirine. The discovery of disubstituted quinoline indicated that the pyridinylamino substituent at the 2-position of quinoline plays an important role in its inhibitory activity against HIV-1 RT. The highly potent HIV-1 RT inhibitors, namely, 4-(2',6'-dimethyl-4'-formylphenoxy)-2-(5″-cyanopyridin-2″ylamino)quinoline (6b) and 4-(2',6'-dimethyl-4'-cyanophenoxy)-2-(5″-cyanopyridin-2″ylamino)quinoline (6d) exhibited half-maximal inhibitory concentrations (IC50) of 1.93 and 1.22 µM, respectively, which are similar to that of nevirapine (IC50 = 1.05 µM). The molecular docking results for these two compounds showed that both compounds interacted with Lys101, His235, and Pro236 residues through hydrogen bonding and interacted with Tyr188, Trp229, and Tyr318 residues through π-π stacking in HIV-1 RT. Interestingly, 6b was highly cytotoxic against MOLT-3 (acute lymphoblastic leukemia), HeLA (cervical carcinoma), and HL-60 (promyeloblast) cells with IC50 values of 12.7 ± 1.1, 25.7 ± 0.8, and 20.5 ± 2.1 µM, respectively. However, 6b and 6d had very low and no cytotoxicity, respectively, to-ward normal embryonic lung (MRC-5) cells. Therefore, the synthesis and biological evaluation of 2-phenylamino-4-phenoxyquinoline derivatives can serve as an excellent basis for the development of highly effective anti-HIV-1 and anticancer agents in the near future.

Biophysical Characterization of Novel DNA Aptamers against K103N/Y181C Double Mutant HIV-1 Reverse Transcriptase.

The human immunodeficiency virus type-1 Reverse Transcriptase (HIV-1 RT) plays a pivotal role in essential viral replication and is the main target for antiviral therapy. The anti-HIV-1 RT drugs address resistance-associated mutations. This research focused on isolating the potential specific DNA aptamers against K103N/Y181C double mutant HIV-1 RT. Five DNA aptamers showed low IC50 values against both the KY-mutant HIV-1 RT and wildtype (WT) HIV-1 RT. The kinetic binding affinity forms surface plasmon resonance of both KY-mutant and WT HIV-1 RTs in the range of 0.06-2 µM and 0.15-2 µM, respectively. Among these aptamers, the KY44 aptamer was chosen to study the interaction of HIV-1 RTs-DNA aptamer complex by NMR experiments. The NMR results indicate that the aptamer could interact with both WT and KY-mutant HIV-1 RT at the NNRTI drug binding pocket by inducing a chemical shift at methionine residues. Furthermore, KY44 could inhibit pseudo-HIV particle infection in HEK293 cells with nearly 80% inhibition and showed low cytotoxicity on HEK293 cells. These together indicated that the KY44 aptamer could be a potential inhibitor of both WT and KY-mutant HIV-RT.

Origin and evolution of HIV-1 subtype A6.

BACKGROUND: HIV outbreaks in the Former Soviet Union (FSU) countries were characterized by repeated transmission of the HIV variant AFSU, which is now classified as a distinct subtype A sub-subtype called A6. The current study used phylogenetic/phylodynamic and signature mutation analyses to determine likely evolutionary relationship between subtype A6 and other subtype A sub-subtypes. METHODS: For this study, an initial Maximum Likelihood phylogenetic analysis was performed using a total of 553 full-length, publicly available, reverse transcriptase sequences, from A1, A2, A3, A4, A5, and A6 sub-subtypes of subtype A. For phylogenetic clustering and signature mutation analysis, a total of 5961 and 3959 pol and env sequences, respectively, were used. RESULTS: Phylogenetic and signature mutation analysis showed that HIV-1 sub-subtype A6 likely originated from sub-subtype A1 of African origin. A6 and A1 pol and env genes shared several signature mutations that indicate genetic similarity between the two subtypes. For A6, tMRCA dated to 1975, 15 years later than that of A1. CONCLUSION: The current study provides insights into the evolution and diversification of A6 in the backdrop of FSU countries and indicates that A6 in FSU countries evolved from A1 of African origin and is getting bridged outside the FSU region.

Sliding of HIV-1 reverse transcriptase over DNA creates a transient P pocket - targeting P-pocket by fragment screening.

HIV-1 reverse transcriptase (RT) slides over an RNA/DNA or dsDNA substrate while copying the viral RNA to a proviral DNA. We report a crystal structure of RT/dsDNA complex in which RT overstepped the primer 3'-end of a dsDNA substrate and created a transient P-pocket at the priming site. We performed a high-throughput screening of 300 drug-like fragments by X-ray crystallography that identifies two leads that bind the P-pocket, which is composed of structural elements from polymerase active site, primer grip, and template-primer that are resilient to drug-resistance mutations. Analogs of a fragment were synthesized, two of which show noticeable RT inhibition. An engineered RT/DNA aptamer complex could trap the transient P-pocket in solution, and structures of the RT/DNA complex were determined in the presence of an inhibitory fragment. A synthesized analog bound at P-pocket is further analyzed by single-particle cryo-EM. Identification of the P-pocket within HIV RT and the developed structure-based platform provide an opportunity for the design new types of polymerase inhibitors.

Structure-Based Design and Discovery of Pyridyl-Bearing Fused Bicyclic HIV-1 Inhibitors: Synthesis, Biological Characterization, and Molecular Modeling Studies.

Two series of new pyridyl-bearing fused bicyclic analogues designed to target the dual-tolerant regions of the non-nucleoside reverse transcriptase inhibitor (NNRTI)-binding pocket were synthesized and evaluated for their anti-HIV activities. Several compounds, such as 6, 14, 15, 21, 30, and 33, were found to be potent inhibitors against the wild-type (WT) HIV-1 strain or multiple NNRTI-resistant strains at low nanomolar levels. Detailed structure-activity relationships were obtained by utilizing the variation of moieties within the corresponding pharmacophores. In vitro metabolic stability profiles and some drug-like properties of selected compounds were assessed, furnishing the preliminary structure-metabolic stability relationships. Furthermore, molecular modeling studies elucidated the binding modes of compounds 6, 15, 21, and 30 in the binding pocket of WT, E138K, K103N, or Y181C HIV-1 RTs. These promising compounds can be used as lead compounds and warrant further structural optimization to yield more active HIV-1 inhibitors.

Using machine learning and big data to explore the drug resistance landscape in HIV.

Drug resistance mutations (DRMs) appear in HIV under treatment pressure. DRMs are commonly transmitted to naive patients. The standard approach to reveal new DRMs is to test for significant frequency differences of mutations between treated and naive patients. However, we then consider each mutation individually and cannot hope to study interactions between several mutations. Here, we aim to leverage the ever-growing quantity of high-quality sequence data and machine learning methods to study such interactions (i.e. epistasis), as well as try to find new DRMs. We trained classifiers to discriminate between Reverse Transcriptase Inhibitor (RTI)-experienced and RTI-naive samples on a large HIV-1 reverse transcriptase (RT) sequence dataset from the UK (n ≈ 55, 000), using all observed mutations as binary representation features. To assess the robustness of our findings, our classifiers were evaluated on independent data sets, both from the UK and Africa. Important representation features for each classifier were then extracted as potential DRMs. To find novel DRMs, we repeated this process by removing either features or samples associated to known DRMs. When keeping all known resistance signal, we detected sufficiently prevalent known DRMs, thus validating the approach. When removing features corresponding to known DRMs, our classifiers retained some prediction accuracy, and six new mutations significantly associated with resistance were identified. These six mutations have a low genetic barrier, are correlated to known DRMs, and are spatially close to either the RT active site or the regulatory binding pocket. When removing both known DRM features and sequences containing at least one known DRM, our classifiers lose all prediction accuracy. These results likely indicate that all mutations directly conferring resistance have been found, and that our newly discovered DRMs are accessory or compensatory mutations. Moreover, apart from the accessory nature of the relationships we found, we did not find any significant signal of further, more subtle epistasis combining several mutations which individually do not seem to confer any resistance.